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ABSTRACT 

We report the discovery of HAT-P-26b, a transiting extrasolar planet orbiting the moderately bright 
V=11.744 Kl dwarf star GSC 0320-01027, with a period P = 4.234516 ± 0.000015d, transit epoch 
T c = 2455304.65122 ± 0.00035 (BJE@), and transit duration 0.1023 ± O.OOlOd. The host star has a 
mass of 0.82 ± 0.03 M , radius of O^Iq';!)" Rq, effective temperature 5079 ± 88K, and metallicity 
[Fe/H] = —0.04 ± 0.08. The planetary companion has a mass of 0.059 ± 0.007 Mj, and radius of 
0.565to °32 Rj yielding a mean density of 0.40 ± 0.10 gem" 3 . HAT-P-26b is the fourth Neptune-mass 
transiting planet discovered to date. It has a mass that is comparable to those of Neptune and Uranus, 
and slightly smaller than those of the other transiting Super-Neptunes, but a radius that is ~65% 
larger than those of Neptune and Uranus, and also larger than those of the other transiting Super- 
Neptunes. HAT-P-26b is consistent with theoretical models of an irradiated Neptune-mass planet 
with a 10 heavy element core that comprises > 50% of its mass with the remainder contained 
in a significant hydrogen-helium envelope, though the exact composition is uncertain as there are 
significant differences between various theoretical models at the Ncptunc-mass regime. The equatorial 
declination of the star makes it easily accessible to both Northern and Southern ground-based facilities 
for follow-up observations. 

Subject headings: planetary systems — stars: individual (HAT-P-26, GSC 0320-01027) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

Transiting exoplanets (TEPs) are tremendously use- 
ful objects for studying the properties of planets out- 
side of the solar system because their photometric tran- 
sits, combined with precise measurements of the ra- 
dial velocity variations of their host stars, enable de- 
terminations of their masses and radii. Of the > 90 
confirmed TEPs discovered to datcFl. all but five are 
Saturn or Jupiter-size gas giant planets with masses 
above 0.1 Mj. The five TEPs below this limit, includ- 
ing the Super-Earths C oRoT-7b (M = 0.0 1 5 ± 0.003 M T 
R = 0.15 ± 0.008 ILeger et alj [2001 IQueloz et all 
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2009), and GJ 12 14b [M = 0.0206 ± 0.003 1 Mj, R = 
0.239 ± 0.012 ICharbonneau et al.1 I2009T) . and the 
Super-Neptunes GJ 436b (M = 0.078 ± 0.006 Mr, R = 
0.376 ± 0.022 R.r IButler et al.ll200l IGillon et all [20071: 
ISouthworthl [20091) . HA T-P-llb (M = 0.08 1 ± 0.009 Mj, 
R = 0.422 ± 0.014i? T : [Bakos et aL|[2010h . and Kepler- 
4b (M = 0.081 ± 0.014 Mj, R 
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Rj', 

IBorucki et all [20T0I: IKipping fc Bakoil2010aD are likely 
composed primarily of elements heavier than hydrogren 
and helium, and are therefore assumed to be qualitatively 
different from the more massive gas giants. In addition 
to these five TEPs, the cand idate TEP Kepler-9 d has an 
estimated radius of 1.4 R @ (IHolman et al.l 120101). a nd is 
most likely a low-mass planet ( Torres et al.ll2010lL but 
currently lacks a mass determination. 

While the gas giant planets exhibit a wide range of 
radii at fixed mass (and hence a great diversity in their 
physical structure at fixed mass), the low mass TEPs, 
together with the six Solar System planets smaller than 
Saturn, appear to follow a nearly monotonic relation be- 
tween mass and radius. The two Super-Neptunes with 
precise radius measurements (GJ 436b and HAT-P-llb) 
have radii that are similar to one another (to within 15%) 
as well as to Uranus (M = 0.0457 Mj, R = 0.358 i?JB) 
and Neptune (M = 0.0540 Mj, R = 0.346 i?j). While the 
mass and radius of K epler-4b given in the discovery paper 
(|Borucki et al.|[2010t) are similar to those of GJ 436b and 

13 Solar system masses are taken from the IAU WG on NSFA 
report of current best estimates to the 2009 IAU General Assem- 
bly retrieved from http://maia.usno.navy.mil/NSFA/CBE.html 
We adopt equatorial r adii tor tne Solar System planets Ironi 
ISeidelmann et al.l H2007T ). 
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HAT-P-llb, a reanalysis by iKipping fc Bakoi (|2010al) 
finds that the radius may be ~ 40% larger, though with 
a 20% uncertainty, it may still be similar to the other 
Super-Neptunes. The lack of significant scatter in the 
radii among Uranus, Neptune, and the Super-Neptunes 
is perhaps surprising given the vast range of radii per- 
mitted by theoretical structure models for planets in 
thi s mass range. For example, the theoretical models 
by iFortnev et al.l (|2007l ) predict that a 1 Gyr nonirradi- 
ated Neptune-mass planet may have a radius that ranges 
from 0.14 i?j (pure iron composition) to 0.29 Rj (pure 
water ice composition) or 0.86 i?j (pure gas composi- 
tion). These same models also predict that the radii of 
gas-dominated Neptune-mass planets should be far more 
sensitive to stellar irradiation than those of Jupiter-mass 
planets. For example, a 1 Gyr pure hydrogen-helium 
Neptune-mass planet at 0.045 AU has a radius of 1.49 Rj 
compared to 1.16 Rj for a similarly irradiated Jupiter- 
mass planet. 

In this paper we present the discovery of HAT-P-26b, a 
TEP orbiting the relatively bright star GSC 0320-01027 
with a mass similar to that of Neptune, but with a radius 
of 0.57 Rj or 65% larger than that of Neptune. This is 
the 26th TEP discovered by the Hun garian-made Auto- 
mated Telescope Network (HATNet: IBakos et al.l 12004ft 
survey. In operation since 2003, HATNet has now cov- 
ered approximately 14% of the sky, searching for TEPs 
around bright stars (8 < r < 14.5). HATNet operates 
six wide-field instruments: four at the Fred Lawrence 
Whipple Observatory (FLWO) in Arizona, and two on 
the roof of the hangar servicing the Smithsonian Astro- 
physical Observatory's Submillimeter Array, in Hawaii. 

The layout of the paper is as follows. In Section [2] we 
report the detection of the photometric signal and the 
follow-up spectroscopic and photometric observations of 
HAT-P-26. In Section [3] we describe the analysis of the 
data, beginning with the determination of the stellar pa- 
rameters, continuing with a discussion of the methods 
used to rule out nonplanetary, false positive scenarios 
which could mimic the photometric and spectroscopic ob- 
servations, and finishing with a description of our global 
modeling of the photometry and radial velocities. Our 
findings are discussed in Section 

2. OBSERVATIONS 

2.1. Photometric detection 

The transits of HAT-P-26b were detected with the 
HAT-5 and HAT-6 telescopes in Arizona, and with the 
HAT-8 and HAT-9 telescopes in Hawaii. Two regions 
around GSC 0320-01027, corresponding to fields inter- 
nally labeled as 376 and 377, were both observed on a 
nightly basis between 2009 Jan and 2009 Aug, whenever 
weather conditions permitted. For field 376 we gathered 
11,500 exposures of 5 minutes at a 5.5 minute cadence. 
Approximately 1500 of these exposures were rejected by 
our photometric pipeline because they yielded poor pho- 
tometry for a significant number of stars. Each im- 
age contained approximately 17,000 stars down to Sloan 
r ~ 14.5. For the brightest stars in the field, we achieved 
a per- image photometric precision of 4 mmag. For field 
377 we gathered 5200 exposures with the same exposure 
time and cadence; we rejected aproximately 700 expo- 
sures. Each image contained approximately 19,000 stars 



down to Sloan r ~ 14.5. We achieved a similar photo- 
metric precision for the brightest stars in this field. 
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Fig. 1.— Unbinned light curve of HAT-P-26 including all 14,500 
instrumental Sloan r band 5.5 minute cadence measurements ob- 
tained with the HAT-5, HAT-6, HAT-8 and HAT-9 telescopes 
of HATNet (see the text for details), and folded with the pe- 
riod P = 4.2345156 days resulting from the global fit described 
in Section[3]|. The solid li ne s hows the "P1P3" transit model fit 
to the light curve (Section l3.3l l, The lower panel shows a zoomed- 
in view of the transit; the dark filled points show the light curve 
binned in phase using a bin-size of 0.002. 

The calibration of the HATNet frames was carried 
out using standard photometric procedures. The cal- 
ibrated images were then subjec ted to star d e tectio n 
and astrometry, as described in iPal fc Bakosl ([2006). 
Aperture photometry was performed on each image 
at the stellar centroids derived from the Two Micron 
All Sky Survey (2MASS; ISkrutskie etall[2006h catalog 
and the individual astrometric solutions. The result- 
ing light curves were decorrelated (cleaned of trends) 
using the Externa l Parameter Decorrelation (EPD; see 
IBakos et al.l I2010D technique in "constan t" mode and 
the T rend Filtering Algorithm (TEA; see IKovacs et al.1 
2005). The light curves were searched for periodic box- 
shaped signals usin g the Box Least-Squares (BLS; see 
IKovacs et al.l I2002T) method. We detected a significant 
signal in the light curve of GSC 0320-01027 (also known 
as 2MASS 14123753+0403359; a = 14 h 12 m 37.44s, 6 = 
+04°03'36.0"; J2000; V=l 1.744 IDroeee et all l200l . 
with an apparent depth of ~ 4.9 mmag (~ 5.5 mmag 
when using TEA in signal- reconstruction mode), and a 
period of P = 4.2345 days (see Figured]). The drop in 
brightness had a first-to-last-contact duration, relative to 
the total period, of q = 0.0242 ± 0.0002, corresponding 
to a total duration of Pq = 2.455 ± 0.025 hr. 

We removed the transits from the combined 376/377 
light curve and searched for additional transiting objects 
using the BLS method; no significant signals were found 
in the data. We also searched the transit-cleaned light 
curve for periodic variations (e.g. due to stella r rota- 
tion) using the Discrete Fourier Transform (e.g. iKurtzl 
1985), and found no coherent variation with an ampli- 
tude greater than 0.4 mmag. 

2.2. Reconnaissance Spectroscopy 

As is routine in the HATNet project, all candidates 
are subjected to careful scrutiny before investing valu- 
able time on large telescopes. This includes spectro- 
scopic observations at relatively modest facilities to es- 
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tablish whether the transit-like feature in the light curve 
of a candidate might be due to astrophysical phenom- 
ena other than a planet transiting a star. Many of these 
false positives are associated with large radial-velocity 
variations in the star (tens of kms -1 ) that are easily rec- 
ognized. 

To carry out this reconnaissance spectroscopy, we 
made use of the Tillinghast Reflector Echelle Spec- 
trograph (TRES; IFtlres A 120081) on the 1.5 m Tilling- 
hast Reflector at FLWO. This instrument provides high- 
resolution spectra which, with even modest signal-to- 
noise (S/N) ratios, are suitable for deriving RVs with 
moderate precision (< 0.3kms _1 ) for slowly rotating 
stars. We also use these spectra to estimate effective 
temperatures, surface gravities, and projected rotational 
velocities of the target star. Using the medium fiber on 
TRES, we obtained two spectra of HAT-P-26 on the 
nights of 2009 Dec 26 and 2009 Dec 27. The spec- 
tra have a resolution of A/AA « 44,000 and a wave- 
length coverage of ~ 3900-8900 A. The spectra were ex- 
tracted and analyzed accordin g to the procedure out - 
lined bv lBuchhave et all (j2010l) and lQuinn et all (|2010t ). 
The individual velocity measurements of 14.62 km s _1 
and 14.81 km s _1 were consistent with no detectable RV 
variation within the measurement precision. Both spec- 
tra were single- lined, i.e., there is no evidence for addi- 
tional stars in the system. The atmospheric parameters 
we infer from these observations are the following: effec- 
tive temperature T c g+ = 5125 ± 125 K, surface gravity 
logg* = 4.5 ± 0.5 (log cgs), and projected rotational ve- 
locity usini = 1.0 ± 1.0 kms . The effective tempera- 
ture corresponds to an early K dwarf. The mean helio- 
centric RV of HAT-P-26 is 7 RV = +14.72 ± 0.10 kms" 1 . 

2.3. High resolution, high S/N spectroscopy 

Given the significant transit detection by HATNet, and 
the encouraging TRES results that rule out obvious false 
positives, we proceeded with the follow-up of this candi- 
date by obtaining high-resolution, high-S/N spectra to 
characterize the RV variations, and to refine the deter- 
mination of the stel lar parameters. F or this we used the 
HIRES instrument (|Vogt et al.lll994fl on the Keck I tele- 
scope located on Mauna Kea, Hawaii, between 2009 Dec 
and 2010 June. The width of the spectrometer slit was 
0'.'86, resulting in a resolving power of A/AA » 55,000, 
with a wavelength coverag e of -3800-8000 A. 

We obtained 12 exposures through an iodine gas ab- 
sorption cell, which was used to superimpose a dense 
forest of I2 lines on the stellar sp ectrum and establish a n 
accurate wavelength fiducial (see IMarcv fc Butlen H992) . 
An additional exposure was taken without the iodine cell, 
for use as a template in the reductions. Relative RVs in 
the solar system baryc e ntric frame were derived as de- 
scribed bv lButler et al.l (|1996fl . incorporating full model- 
ing of the spatial and temporal variations of the instru- 
mental profile. The RV measurements and their uncer- 
tainties are listed in Table [0 The period-folded data, 
along with our best fit described below in Section 02 are 
displayed in Figure [2j 

In the same figure we show also the Shk index, which 
is a measure of the chromospheric activity of the star 
derived from the flux in the cores of the Ca II H 
and K lines (Figure [3] shows representative Keck spec- 
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Fig. 2. — Top panel: Kcck/HIRES RV measurements for 
HAT-P-26 shown as a function of orbital phase, along with our 
best-fit model (see Table . Zero phase corresponds to the time 
of mid-transit. The center-of-mass velocity has been subtracted. 
Second panel: Velocity O — C residuals from the best fit. The 
error bars include a component from astrophysical/instrumental 
jitter ( 1.6m s" 1 ) added in quadrature to the formal errors (see 
Section l3.3l l. Third panel: Bisector spans (BS), with the mean 
value subtracted. The mea surement from the template spectrum 
is included (see Section [32} . Bottom panel: Chromospheric activ- 
ity index Shk measured from the Keck spectra. Note the different 
vertical scales of the panels. Observations shown twice are repre- 
sented with open symbols. 
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Fig. 3. — Calcium K (left) and H (right) line profile in selected 
Keck/HIRES observations of HAT-P-26. Both panels show three 
spectra overlaid; data taken at high, median, and low activity, as 
characterized by the Sjjk index. Low-level emission is clearly de- 
tected in the cores of both lines. The lack of apparent variation 
(the three plotted spectra are indistinguishable) indicates the chro- 
mospheric stability of this star over the course of our observations. 
The spectra are matched to a common flux/ wavelength scale using 
points outside the H and K line cores. The vertical axes give the 
counts in units of CCD e~ per wavelength bin for the reference 
spectrum. 

tra including the H and K lines for HAT-P-26). This 
index was computed and calibrated to the scale of 
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iVaughan. Preston fc Wilson! (11978D fol l owing the proce- 
dure described by llsaacson fc Fischerl (gOlfl). We find 
a median value of Shk = 0.182 with a standard devi- 
ation of 0.004. Assuming B — V = 0.89 based on the 
effective temperature measure d in Section 13. 1| t his cor- 
responds to log R^y = -4 .992 (jNoves et al.lll984h . From 
llsaacson fc Fischerl (|2010D the lower tenth percentile Shk 
value among California Planet Search (CPS) targets with 
B — V = 0.89 is 0.168. The measured Shk value is only 
slightly higher than this, implying that HAT-P-26 is a 
chromospherically quiet star. We do not detect any sig- 
nificant variation of the Shk index correlated with or- 
bital phase; such a correlation might have indicated that 
the RV variations could be due to stellar activity, casting 
doubt on the planetary nature of the candidate. 

2.4. Photometric follow-up observations 
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Fig. 4. — Unbinned instrumental Sloan i band light curves, ac- 
quired with KeplerCam at the FLWO 1.2 m telescope. The light 
curves have been EPD- and TFA-processed, as described in S 13.31 
The dates of the events are indicated. Curves after the first are 
displaced vertically for c larit y. Our best fit from the global model- 
ing described in Section 13.31 is shown by the solid lines. Residuals 
from the fits are displayed at the bottom, in the same order as the 
top curves. The error bars represent the photon and background 
shot noise, plus the readout noise. 

In order to permit a more accurate modeling of the 
light curve, we conducted additional photometric ob- 
servations with the KeplerCam CCD camera on the 
FLWO 1.2 m telescope. We observed five transit events 
of HAT-P-26 on the nights of 2010 Jan 05, Mar 31, Apr 
04, May 08 and May 25 (Figure H}. These observations 
are summarized in Table [5J 

The reduction of these images, including basic cali- 
bration, astrometry, and aperture photometry, was per- 
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Fig. 5. — The combined Sloan i band KeplerCam light curve 
of HAT-P-26 folded with the period P = 4.2345156 days resulting 
from the global fit described in Section [3] The dark filled cir- 
cles show the light curve binned in folded time with a bin-size 
of 0.005 days. The median uncertainty on the binned points is 
0.19 mmag. 



formed as described by iBakos et alj (|2010f ). We per- 
formed EPD and TFA to remove trends simultaneously 
with the light curve modeling (for more details, see 
Section [31 and IBakos et al.l I2010T ) . The final time series 
are shown in the top portion of Figure 21 along with our 
best-fit transit light curve model described below; the 
individual measurements are reported in Table [3l The 
combined phase-folded follow-up light curve is displayed 
in Figure [5] 

3. ANALYSIS 

3.1. Properties of the parent star 

Fundamental parameters of the host star HAT-P-26 
such as the mass (M*) and radius (-R*), which are needed 
to infer the planetary properties, depend strongly on 
other stellar quantities that can be derived spectroscop- 
ically. For this we have relied on our template spec- 
trum obtained with the Keck/HIRES instrument, and 
the an alysis package known as S pectroscopy Made Easy 
(SME; iValenti fc Piskunov 1996), a long w ith the atomic 
line database of IValenti fc Fischerl (|2005[ h SME yielded 
the following values and uncertainties: effective tem- 
perature Teg* = 5079 ± 88 K, stellar surface gravity 
log#* = 4.53 ± 0.06 (cgs), metallicity [Fe/H] = -0.04 ± 
0.08 dex, and projected rotational velocity v sinz = 1.8 ± 
0.5 km s _1 , in which the uncertainties for T e ff+ and [Fe/H] 
have been increased by a factor of two over their formal 
values to include our estimates of the systematic uncer- 
tainties. 

In principle the effective temperature and metallicity, 
along with the surface gravity taken as a luminosity in- 
dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. However, the effect of log g+ on the spectral 
line shapes is quite subtle, and as a result it is typ- 
ically difficult to determine accurately, so that it is a 
rather poor luminosity indicator in practice. For plan- 
etary transits a stronger constraint is often provided 
by p* the mean stellar density, which is closely related 
to the a/R+ normalized semimajor axis. The quan- 
tity p+ can be derived d irectly from the combination of 
the transit light curves (jSeager fc Mallen-Ornelai f2003; 
ISozzetti et al.ll2007f ) and the RV data (required for eccen- 



tric cases, see Section 1373)) . This, in turn, allows us to im- 
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TABLE 1 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-26. 



BJD 


RV a 


<JRV b 


BS 


CBS 


Shk c 


Phase 


(2,454,000+) 


(ms- 1 ) 


(ms- 1 ) 


(ms" 1 ) 


(m s -1 ) 






1193.11925 


6.67 


1.87 


3.71 


2.51 


0.1870 


0.661 


1193.12855 






4.74 


3.90 


0.1830 


0.663 


1194.16006 


10.75 


1.62 


-1.85 


3.76 


0.1850 


0.907 


1252.02017 


2.23 


1.89 


8.33 


3.95 


0.1890 


0.571 


1285.14629 


-7.14 


2.58 


-24.47 


4.35 


0.1760 


0.394 


1320.84584 


8.96 


2.00 


-13.90 


5.02 


0.1840 


0.824 


1343.78649 


-3.54 


1.79 


-8.59 


3.23 


0.1870 


0.242 


1350.92272 


4.54 


2.18 


-14.55 


4.65 


0.1740 


0.927 


1351.91410 


-6.17 


1.70 


7.80 


4.04 


0.1800 


0.161 


1372.77275 


-6.44 


1.65 


4.63 


2.63 


0.1800 


0.087 


1373.75770 


-9.14 


1.71 


5.32 


4.79 


0.1810 


0.320 


1374.91646 


2.06 


1.78 


22.23 


4.36 


0.1780 


0.593 


1375.80851 


6.63 


1.71 


6.59 


3.90 


0.1820 


0.804 



Note. — For the iodine-free template exposures there is no RV mea- 
surement, but the BS and Sjjk index can still be determined. 
a The zero-point of these vel ociti es is arbitrary. An overall offset 7 re i fitted 
to these velocities in Section l3.3l has not been subtracted. 
b Internal errors excludin g th e component of astrophysical/instrumental 
jitter considered in Section 13.31 

c "Shk chromospheric activity index, calibrated to the scale of 
Vaughan, Preston & Wilson] (fl97a) . 



TABLE 2 

Summary of photometric follow-up observations 



Facility 



Date 



Number of Images Median Cadence (s) 



Filter 



KeplerCam/FLWO 1.2 m 
KeplerCam/FLWO 1.2 m 
KeplerCam/FLWO 1.2 m 
KeplerCam/FLWO 1.2 m 
KeplerCam/FLWO 1.2 m 



2010 Jan 05 
2010 Mar 31 
2010 Apr 04 
2010 May 08 
2010 May 25 



191 
161 
291 

596 
298 



40 


Sloan i 


band 


59 


Sloan i 


band 


61 


Sloan i 


band 


11 


Sloan i 


band 


59 


Sloan i 


band 



TABLE 3 

High-precision differential photometry of HAT-P-26 



BJD 


Mag a 


o-Mag Mag(orig) b Filter 


(2,400,000+) 






55202.94802 


-0.00279 


0.00105 10.17790 i 


55202.94847 


0.00137 


0.00106 10.18150 i 


55202.94909 


0.00132 


0.00106 10.18130 i 


55202.94954 


0.00079 


0.00106 10.18110 i 


55202.95018 


-0.00047 


0.00105 10.17920 i 


55202.95063 


-0.00029 


0.00106 10.17980 i 


55202.95126 


-0.00131 


0.00106 10.17960 i 


55202.95171 


-0.00097 


0.00105 10.17950 i 


55202.95237 


0.00334 


0.00107 10.18530 i 


55202.95282 


-0.00366 


0.00106 10.17800 i 



Note. — This table is available in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form 
and content. 



a The out-of-transit level has been subtracted. These magnitudes have 
been subjected to the EPD and TFA procedures, carried out simulta- 
neously with the transit fit. 

15 Raw magnitude values without application of the EPD and TFA 
procedures. 

prove on the determination of the spectroscopic param- 
eters by supplying an indirect constraint on the weakly 
determined spectroscopic value of log that removes de- 
generacies. We take this approach here, as described be- 
low. The validity of our assumption, namely that the 
adequate physical model describing our data is a plan- 
etary transit (as opposed to a blend), is shown later in 



Section l3~2l 

Our values of T e ff*, logg*, and [Fe/H] were used to 
determine auxiliary quantities needed in the global mod- 
eling of the follow-up photometry and radial velocities 
(specifically, the limb-darkening coefficients). This mod- 
eling, the details of which are described in Section l3~3l 
uses a Monte Carlo approach to deliver the numerical 
probability distribution of p* and other fit ted variables . 
For further details we refer the reader to iPall ((2009b) . 
When combining (a luminosity proxy) with assumed 
Gaussian distributions for T e ff* and [Fe/H] based on the 
SME determinations, a comparison with stellar evolu- 
tion models allows the probability distributions of other 
stellar properties to be inferred, including log <?*. Here 
we use the stellar ev olution calculat ions from the Yonsei- 
Yale (YY) series bv IYi et all IpOOl . 

For the case of HAT-P-26b, the eccentricity is poorly 
constrained by the RV data due to the low semiamplitude 
of the signal. This in turn leads to a significant uncer- 
tainty on p+. However, not all combinations of [Fe/H], 
T G ff* and p± are realized by physical stellar models. In 
particular, if we conservatively adopt a maximum stel- 
lar age of 13.8 Gyr, c orresponding approxi mately to the 
age of the universe (jKomatsu et alJ[2CM find 13.75 ± 
0.11 Gyr), and a minimum age of 100 Myr, correspond- 
ing roughly to the zero-age main-sequenc43, stars with 

14 The lack of evidence for stellar activity indicates that HAT- 
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T eff * = 5079 and [Fe/H] = -0.04 are not found to have 
densities in the range 0.24 gem" 3 < < 2.06 gem" 3 
or surface gravities in the range 3.915 < log 5* < 4.514 
(here 0.24 gem" 3 corresponds to an evolved star with 
M = 0.94 Mq, while 2.06 gem" 3 corresponds to a main 
sequence star with M = 0.79 M Q ). Figure|5]shows the in- 
ferred location of the star in a diagram of p+ versus T e ff*, 
analoguous to the classical H-R diagram, for three cases: 
fixing the eccentricity of the orbit to zero, allowing the 
eccentricity to vary, and allowing the eccentricity to vary, 
but only accepting parameter combinations which match 
to a position in the YY isochrones with 0.1 Gyr < age < 
13.8 Gyr. The stellar properties and their approximate 
la and 2a con fidence boundar ies are displayed against 
the backdrop of lYi et all ()2001[ ) isochrones for the mea- 
sured metallicity of [Fe/H] = —0.04, and a range of ages. 
For the zero eccentricity case the comparison against 
the model isochrones was carried out for each of 30,000 
Monte Carlo trial sets (see Section 153)) . We find good 
overlap between the trials and the model isochrones-in 
71% of the trials, the [Fe/H], T e s* and P* parameter com- 
bination matched to a physical location in the H-R dia- 
gram that has an age that is within the aforementioned 
range. However, when the eccentricity is allowed to vary, 
the model for the light curves and RV data tends toward 
low values of which may only be fit by pre-main se- 
quence stellar models, or stellar models older than the 
age of the universe. In this case only 15% out of 100,000 
Monte Carlo trial sets match to physical locations in the 
H-R diagram with ages within the allowed range. By 
requiring the star to have an age between 0.1 Gyr and 
13.8 Gyr, we effectively impose a tighter constraint on 
the orbital eccentricity than is possible from the RV data 
alone (we find an eccentricity of e = 0.124±0.060, as com- 
pared with e = 0.24 ± 0.12 when not requiring a match 
to the stellar models; see also Section l3~3| . 

Adopting the parameters which result from allowing 
the eccentricity to vary while requiring the star to have 
an age between 0.1 Gyr and 13.8 Gyr yields a stellar sur- 
face gravity of log g+ = 4.56 ±0.06, which is very close to 
the value from our SME analysis. The values for the at- 
mospheric parameters of the star are collected in TablelH 
together with the adopted values for the macroturbulent 
and microturbulent velocities. 

The stellar evolution modeling provides color indices 
that may be compared against the measured values as 
a sanity check. The best available measurements are 
the ne ar-infrared magnitud es from the 2MASS Cata- 
logue (jSkrutskie et al.l l2006f). J 2 mass = 10.080 ± 0.022, 
#2MASS = 9.685 ± 0.023 and K 2M ass = 9.581 ± 
0.023; which we have converted to the photometric sys- 
tem of the mode ls (ESO) using the transformations by 
ICarpenterl (|2001[ ) . The resulting measured color index is 
J-K = 0.530 ±0.035. This is within la of the predicted 
value from the isochrones of J — K = 0.55 ± 0.02. The 
distance to the object may be computed from the abso- 
lute K magnitude from the models (Mr = 3.98 ± 0.19) 
and the 2MASS K s magnitude, which has the advantage 
of being less affected by extinction than optical magni- 
tudes. The result is 134lg 8 pc, where the uncertainty ex- 
cludes possible systematics in the model isochrones that 

P-26 is unlikely to be a pre-main sequence star. 
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Fig. 6. — Model isochrones from lYi et al.l 11200 If ) for the mea- 
sured metallicity of HAT-P-26, [Fe/H]= -0.04, and ages in 2 Gyr 
increments between 1 and 13 Gyr (left to right). Note that a log- 
arithmic scale is used for the vertical axis. The measured values 
of Teff* and p* are shown together with their approximate la and 
2<r confidence boundaries for models where the eccentricity is fixed 
to (top), the eccentricity is allowed to vary (bottom), and the 
eccentricity is allowed to vary, but only models which match to an 
isochrone with an age between 0.1 Gyr and 13.8 Gyr are accepted 
(bottom). We adopt the parameters shown in the bottom panel. 

are difficult to quantify. 

An additional check on our stellar model can be per- 
formed by comparing the isochrone-based age estimate to 
activity-based age estimates. U sing the activity-rotation 
and a ctivity- age relations from [Mama i ok "&: Hillenbrand! 
(2008, equations 5 and 3) w e c onvert the value of 
log i?H K determined in Section 12.31 into a Rossby num- 
ber (Ro = Prot/i~ c , where P ro t is the rotation period 
and r c is the convective turnover time-scale) and an age. 
We find R = 2.2 ± 0.2, and log(r) = 9.80 ± 0.15 or 
t = 6.4+^ 9 Gyr, where we adopt the estimated uncer- 
tainties on Ro and log(r) from lMamaiek & Hillen brand! 
p008l) . The R value may be converted to a rotation 
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TABLE 4 

Stellar parameters for HAT-P-26 



Parameter Value Source 



Spectroscopic properties 

T eff * (K) 5079 ± 88 SME a 

[Fe/H] -0.04 ± 0.08 SME 

vsini (kms" 1 ) . . . 1.8 ± 0.5 SME 

Umac (kms" 1 ).... 2.95 SME 

v mic (kms" 1 ) .... 0.85 SME 

7rv (kms" 1 ) +14.72±0.10 TRES 

Photometric properties 

V (mag) 11.744 TASS 

V-I c (mag) 0.96 ±0.11 TASS 

J (mag) 10.080 ± 0.022 2MASS 

H (mag) 9.685 ± 0.023 2MASS 

K s (mag) 9.581 ± 0.023 2MASS 

Derived properties 

M 4 (Mq) 0.816 ±0.033 YY+a/i?*+SME b 

R+ (Rq) 0.788to°tl YY+a/R*+SME 

logg* (cgs) 4.56 ±0.06 YY+a/R*±SME 

L* (L ) 0.381°'^ YY+a/i?*±SME 

(mag) 6.03 ±0.24 YY+a/2t.*+SME 

M K (mag,ESO) . . 3.98 ±0.19 YY+a/i?*±SME 

Age (Gyr) S.oj^g YY+a/i?*+SME 

Distance (pc) 134;tg 8 YY+a/i?*±SME 



SME — "Spectroscopy Made Easy" package for the analysis of 



high-resolution spectra (Valcnti & Piskunov 1996). These parame- 
ters rely primarily on SME, but have a small dependence also on 
the iterative analysis incorporating the isochronc search and global 
modeling of the data, as described in the text. 

b YY+a/fl,+SME = Based on the YY isochrones lIYi et al.H200in . 
a/Ft* as a luminosity indicator, and the SME results. 

period using the relation for r c given by iNoves et al.1 
(1984). We find P rot = 48 ± 4d. The rotation pe- 
riod and color may also be used to obtain a separate 
age estimate from the gyrochro nology relation given by 
iMamaiek fe Hillenbrand! (|2008l equation 12). This gives 
t = 7.8^^2 Gyr, where the scatter on this relation is less 
than the scatter on the age-activity relation because it 
includes a correction for stellar color. The age inferred 
from log i?HK 1S consistent with the isochrone-based age 
of 9.0i"4'g Gyr. The equatorial rotation period inferred 
from the spectroscopic determination of v sin i assuming 
sini = 1 is P ro t, vsini = 22.3^4 7 6 2 d, which is shorter than 
the expected value based on logi?^ K , though the upper 
limit is poorly constrained. 

As discussed below in Section [3.31 we measure 
a RV jitter of 1.6 ms " 1 for HAT-P-26. From 
llsaacson fc Fischer] (|2010f ) the lower tenth percentile jit- 
ter among HIRES/Keck observations for CPS stars with 
0.7 < B-V < 1.0 and Shk = 0.182 is 2.17ms" 1 , imply- 
ing that HAT-P-26 has an exceptionally low jitter value- 
only a handful of stars in this color range have measured 
jitter values less than that of HAT-P-26. We note that 
the jitter may be higher (2.4 m s" 1 ) if the orbit is circular, 
though this value is still quite low. 

3.2. Excluding blend scenarios 

Our initial spectroscopic analyses discussed in 
Section |2"T21 and Section |2~3"1 rule out the most obvious as- 
trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 



resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric and spectroscopic signatures we see. In the fol- 
lowing sections we consider and rule out the possibility 
that such scenarios may have caused the observed pho- 
tometric and spectroscopic features. 

3.2.1. Spectral line-bisector analysis 

Following lQueloz et all ()2001[) ; lTorres et all (|2007l ). we 
explored the possibility that the measured radial veloc- 
ities are not real, but are instead caused by distortions 
in the spectral line profiles due to contamination from a 
nearby unresolved eclipsing binary. A bisector analysis 
b ased on the Keck sp ectra was done as described in §5 
of Ba kos et al.l ([2007D . While the bisector spans show no 
evidence for variation in phase with the orbital period, 
the scatter on these values (13.3 ms" 1 RMS) is large rel- 
ative to the RV semi- amplitude (~ 8 ms" 1 ), and thus the 
lack of variation does not provide a strong constraint on 
possible blend scenarios. We note that some of the scat- 
ter in the bisector spans may be due to contamination 
from the sky background (predominately moonlight- 
correcting t he bisector span s for s ky contamination as 
discussed by lHartman et all (|2010D reduces the RMS to 
~ 10.0ms" 1 , however the precision is still insufficient to 
rule out blend scenarios. 

3.2.2. Contamination from a background eclipsing binary 

Following o ur earlier work ()Bakos et al.l 120101 : 

lHartman et al.l I2009D we make use of the high proper 
motion of HAT-P-26 to rule out the possibility that the 
observed transits and RV variation may be due to a 
background eclipsing binary that is aligned, by chance, 
with the foreground Kl dwarf HAT-P-26. To reproduce 
the observed ~ 0.6% deep transit, the background object 
cannot be more than 5.6 mag fainter than HAT-P-26 
(objects fainter than this would contribute less than 
0.6% of the total combined light and so could not cause 
the transit even if they were to be completely eclipsed by 
an object that emits no light). Because HAT-P-26 has a 
high proper motion (148.5 ± 2.7masyr~ 1 ; iRoeser et al.1 
120101) it is possible to use the Palomar Observatory Sky 
Survey plates from 1950 (POSS-I, red and blue plates) 
to view the sky at the current position of HAT-P-26. 
Between 1950 and the follow-up observations in 2010, 
HAT-P-26 has moved ~ 8'/9. Figure [7] shows an image 
stamp from the POSS-I plate compared with a recent 
observation from the FLWO 1.2 m. We can rule out a 
background object down to R ~ 19 mag within ~ 3" 
of the current position of HAT-P-26. Any background 
object must be > 7.5 mag fainter than HAT-P-26 and 
thus could not be responsible for the observed transit. 

3.2.3. Detailed blend modeling of a hierarchical triple 

F ollowing iBakos et all (|2010l ). lHartman etldl (|2009D . 
and iTorres et al.l (|2005l) we attempt to model the obser- 
vations as a hierarchical triple system. We consider 4 
possibilities: 

1. One star orbited by a planet, 

2. Three stars, 2 fainter stars are eclipsing, 

3. Two stars, 1 planet, planet orbits the fainter star, 
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Fig. 7. — Images of a 2' X 2' field containing HAT-P-26 from the POSS-I Red survey (left), and from our FLWO 1.2m i band follow-up 
observations (right, see Section l2,4jl . North is up and east is to the left in both images. The dates of the exposures are 1950 April 21 
and 2010 May 08, respectively. The cross marks the position of HAT-P-26 in 1950 and the triangle marks the position in 2010. Between 
these two dates HAT-P-26 has moved ~ 8'.'9 to the southeast. From the POSS-I image, we can rule out the presence of stars brighter than 
R ~ 19 at the current position of HAT-P-26. 
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Fig. 8.— Left: \ 2 for the scenario of a binary star with a planet orbiting the fainter star (case 3 in Section 13. 2. 3|l vs. the mass of the 
planet host. Points above the dotted horizontal line are rejected at the 3er confidence level. The mass of the brightest star in the system is 
fixed to 0.788 Mq from the spectroscopically determined effective temper ature, surface gravity and metallicity. Right: \ 2 f° r the scenario 
of a binary star with a planet orbiting the brighter star (case 4 in Section |3. 2. 3D vs. the mass of the fainter star. We also show the radius 
of the planet (right axis). Points on the \ 2 curve above the dotted horizontal line are rejected at the 3er confidence level. 



4. Two stars, 1 planet, planet orbits the brighter star. 

Here case 1 is the fiducial model to which we compare the 
various blend models. We model the observed follow-up 
and HATNet light curves (including only points that are 
within one transit duration of the primary transit or sec- 
ondary eclipse assuming zero eccentricity) together with 
the 2MASS and TASS photometry. In each case we fix 
the mass of the brightest star to 0.788 M©; this ensures 
that we reproduce the effective temperature, metallic- 
ity, and surface gravity determined from the SME anal- 
ysis when using the Padova isochrones (see below). We 
have also attempted to perform the fits described be- 
low allowing the mass of the brightest star to vary. We 
find that in this case the mass of the brightest star is 
still constrained to be close to 0.788 M by the broad- 
band photometry, even if the spectroscopic parameters 
are not included. We therefore conclude that fixing the 



mass of the brightest star is justified. In all cases we 
vary the distance to the system and two parameters al- 
lowing for dilution in the two HATNet light curves, and 
we include simultaneous EPD and TFA in fitting the 
light curves (see Section 1373]) . In each case we draw the 
stellar rad ii and magnitudes from the 13.0 Gyr Padova 
isochronc (Gir ardi et al.ll2000D . extended below 0.15 Mq 
with the iBaraffe et al.l (|1998l) isochrones. We use these 
rather than the YY isochrones for this analysis because 
of the need to allow for stars with M < 0.4 Mq, which is 
the lower limit availa ble for the YY models. W e use the 
JKTEBOP program (jSouthworth et al.l l2004al lbl) which 
is based on the Eclips i ng Binary Orbit Program (EBOP ; 
IPopper fc E"tiel[l98l lEtiel[T98lt iNelson fc Davisl[l972T) 
to generate the model light curves. We optimize the 
free parameters using the Downhill Simplex Algorithm 
together with classical linear least squares for the EPD 
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and TFA parameters. We rescale the errors for each light 
curve such that \ 2 P er degree of freedom is 1.0 for the 
out of transit portion of the light curve. Note that this 
is done prior to applying the EPD/TFA corrections, as a 
result x 2 P er degree of freedom is less than 1.0 for each 
of the best-fit models discussed below. If the rescaling is 
not performed, the difference in x 2 between the best-fit 
models is even more significant than what is given be- 
low, and the blend-models may be rejected with higher 
confidence. 

Case 1: 1 star, 1 planet: In addition to the parame- 
ters mentioned above, in this case we vary the radius of 
the planet and the impact parameter of the transit. The 
best-fit model has Xcasci = 3140.1 for 3567 degrees of 
freedom. The parameters that we obtain for the planet 
are comparable to those obtained from the global mod- 
elling described in Section 13.31 

Case 2: 3 stars: For case 2 we vary the masses of 
the eclipsing components, and the impact parameter of 
the eclipse. We find no model of three stars which re- 
produces the observations. The transit depth and dura- 
tion cannot be fit when the three stars are constrained 
to fall on the same isochrone, and the brightest star 
has M = 0.788 Mq. The best-fit case 2 model con- 
sists of equal masses for the bright est two stars, and 

0. 08 Mq (the lowest stellar mass in the Baraffc et al.ll998l 
isochrones) for the transiting star. Such a model is in- 
consistent with our spectroscopic observations (it would 
have been easily identified as a double-lined binary at one 
of the quadrature phases), and as we will show, can be 
rejected from the light curves alone. The best-fit case 2 
model yields Xc aS e2 = 3310.2 for 3566 degrees of freedom 
and produces model transits that are too deep compared 
to the observed transit. The case 1 model achieves a 
lower x 2 with fewer parameters than the case 2 model, 
so the case 1 model is preferred over the case 2 model. 
Assuming that the errors are uncorrelated and follow a 
Gaussian distribution, the case 2 model can be rejected in 
favor of the case 1 model at the > 116cr confidence level. 
Alternatively, one might suppose that any apparent cor- 
relations in the residuals of the best-fit case 2 model are 
not due to errors in the model but instead are due to un- 
corrected systematic errors in the measurements; large 
systematic errors in the measurements would increase the 
probability that case 1 might give a better fit to the data, 
by chance, than case 2. To establish the statistical signf- 
icance with which we may reject case 2 while allowing for 
possible systematic errors in the measurements, we con- 
duct 1000 Monte Carlo simulations in which we assume 
the best-fit case 2 model scenario is correct, shuffle the 
residuals from this fit in a manner that preserves the cor- 
relations (this is done by taking the Fourier Transform 
of the residuals, randomly changing phases of the Trans- 
form while preserving the amplitudes, and transforming 
back to the time domain), and fit both the case 2 and 
case 1 models to the simulated data. The median value of 
Xcasc2 — Xcasci ^ s — 17.6 with a standard deviation of 26.3. 
None of the 1000 trials have Xc asC 2 ~ Xcasci > 170.1, the 
measured value. Based on this analysis we reject case 2, 

1. e. the hierarchical triple star system scenario, at the 7a 
level. 

Case 3: 2 stars, planet orbits the fainter star. In this 
scenario HAT-P-26b is a transiting planet, but it would 



have a radius that is larger than what we infer (it may be 
a Saturn- or Jupiter-size planet rather than a Neptune- 
size planet). For this case we vary the mass of the faint 
planet-hosting star, the radius of the planet, and the im- 
pact parameter of the transit. We assume the mass of 
the planet is negligible relative to the mass of its faint 
host star. Figure [S] shows x 2 as a function of the mass 
of the planet-hosting star for this scenario. The best-fit 
case 3 solution has x 2 = 3147.6 for 3566 degrees of free- 
dom, and corresponds to a system where the two stars 
are of equal mass and the planet has a radius of 0.8 
As the mass of the planet host is decreased the value of 
X 2 increases. Repeating the procedure outlined above to 
establish the statistical significance at which we may re- 
ject case 3 we find that the 3cr limit on Xcn SC 3 ~ Xcasci ^ s 
9.0, which results in a 3cr lower-limit of 0.77 Mq on the 
mass of the planet hosting star. We may thereby place 
a 3(7 lower limit of 0.74 on the y-band luminosity ratio 
of the two stars. A second set of lines with a luminos- 
ity ratio of > 0.74 would have easily been detectable in 
both the Keck and TRES spectra unless the stars had 
very similar 7 velocities (the spectral lines are quite nar- 
row with vsini = 1.8 ± 0.5kms _1 ). The poor fit for 
this blend model relative to the fiducial model together 
with the tight constraints on the relative 7 velocities and 
luminosity ratios of the stars in the blend models that 
may yet fit the data leads us to reject this blend scenario 
in favor of the simpler model of a single star hosting a 
transiting planet. 

Case 4 : % stars, planet orbits the brighter star. As in 
the previous case, in this scenario HAT-P-26b is a tran- 
siting planet, but the dilution from the blending star 
means that the true radius is larger than what we infer 
in Case 1. In this case we vary the radius of the planet, 
the mass of the faint star, and the impact parameter of 
the transit. Again we assume the mass of the planet 
is negligible relative to the mass of its bright host star. 
Figure [8] shows x 2 as a function of the mass of the faint 
star. The smallest value of % 2 is achieved when the faint 
star contributes negligible light to the system, which is 
effectively equivalent to the fiducial scenario represented 
by Case 1. Two effects cause % 2 to increase with stellar 
mass. First, the shape of the transit is subtly changed 
in a manner that gives a poorer fit to the observations. 
Second, when the mass of the faint star is less than that 
of the transit host the model broad-band photometry for 
the blended system is redder than for the single star sce- 
nario, and is inconsistent with the observed photometry. 
This gives rise to the peak in % 2 at M ~ 0.65 Mq. The 
case 4 model where the faint star has M > 0.77 Mq can 
be rejected as in Case 3. For lower masses, we place a 
3(7 upper limit of 0.55 Mq on the mass of the faint star, 
yielding a 3<r upper limit on the luminosity ratio of 0.1. 
We conclude that at most the planet radius R p may be 
8% larger than what we find in Section |3~31 if there is an 
undetected faint secondary star in the system. 

3.3. Global modeling of the data 

Here we summarize the procedure that we followed 
to model the HATNet photometry, the follow-up pho- 
tometry, and the radial velocities simultan eously. This 
proced ure is described in greater detail in iBakos et al.l 
(2010). The follow-up ligh t curves were modeled using 
analytic formulae based on lMandel fc Agol f2002), with 
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quadratic limb darkening coefficie nts for the Sloa n i band 
interpolated from the tables by iClaretl (|2004j ) for the 
spectroscopic parameters of the star as determined from 
the SME analysis (Section EH]). We mo deled the HAT- 
Net d ata using an approximati on to the Mandel & Agol 
( 2002) formulae as described in IBakos et al l (|2010[) . The 
RVs were fitted with an eccentric Keplerian model. 

Our physical model consisted of 8 main parameters, 
including: the time of the first transit center observed 
with HATNet (taken to be event -74), T c _ 74 , and that 
of the last transit center observed with the FLWO 1.2 m 
telescope, T c + 4o, the normalized planetary radius p = 
Rp/R*, the square of the impact parameter 6 2 , the re- 
ci procal of the h alf du ration of the transit (/R+ as given 
in iBakos et al.l (|2010f) . the RV semiamplitude K : and 
the Lagrangian elements k = ecosw and h = esinw, 
where w is the longitude of periastron. Five additional 
parameters were included that have to do with the in- 
strumental configuration. These are the HATNet blend 
factors -Binst,376, and -Binst,377, which account for possi- 
ble dilution of the transit in the HATNet light curves 
from background stars due to the broad PSF (~ 24" 
FWHM), the out-of-transit magnitudes for each HATNet 
field, M ,HATNct,376 an d M),HATNet,377, and the relative 
zero-point 7 rc i of the Keck RVs. The physical model was 
extended with an instrumental model for the follow-up 
light curves that describes brightness variations caused 
by systematic errors in the measurements. We adopted 
a "lo cal" EPD- and "global" TFA-model (jBakos et al.l 
l20ioT) . using 20 template stars for the TFA procedure 
and six EPD parameters for each follow-up light curve. 
In summary, the total number of fitted parameters was 13 
(physical model with 5 configuration-related parameters) 
+ 30 (local EPD) + 20 (global TFA) = 63, i.e., much 
smaller than the number of data points (1450, counting 
only RV measurements and follow-up photometry mea- 
surements). 

As described in IBakos et al.l ()2010D . we use a combi- 
nation of the dow nhill simplex method (AMOEBA; see 
iPress et al.1 11992T ). the classical linear least squares al- 
gorithm, and the Mark ov Chain Monte-Carlo method 
(MCMC, see iFordl [2006] ) to obtain a best-fit model to- 
gether with a posteriori distributions for the fitted pa- 
rameters. These distributions were then used to obtain 
a posteriori distributions for other quantities of interest, 
such as p+. As described in Section l3.il /?* was used 
together with stellar evolution models to obtain a poste- 
riori distributions for stellar parameters, such as Af* and 
i?*, which are needed to determine M p and R p . 

The resulting parameters pertaining to the light curves 
and velocity curves, together with derived physical pa- 
rameters of the planet, are listed under the "Adopted 
Value" column heading of Table [5j Included in this ta- 
ble is the RV "jitter". This is a component of assumed 
astrophysical noise intrinsic to the star, possibly with a 
contribution from instrumental errors as well, that we 
added in quadrature to the internal errors for the RVs 
in order to achieve x 2 /dof = 1 from the RV data for the 
global fit. Auxiliary parameters not listed in the table 
are: T c _ 74 = 2454860.02709 ± 0.00147 (BJD), T c +40 = 
2455342.76185 ± 0.00041 (BJD), the blending factors 

#instr,376 = °- 92 ± °- 05 > and Anstr.377 = 0.84 ± 0.10 for 

the HATNet field 376 and 377 light curves, respectively, 



and 7 ro i = — 2.8 ± 0.8m s _1 . The latter quantity rep- 
resents an arbitrary offset for the Keck RVs, but does 
not correspond to the true center-of-mass velocity of the 
system, which was listed earlier as 7rv m Tabled 

We find a mass for the planet of M p = 0.059±0.007Mj 
and a radius of R p = 0.565^o;q32 Rj, leading to a mean 
density p p = 0.40 ± 0.10 gem -3 . We also find that the 
eccentricity of the orbit may be different from zero: e = 
0.124 ± 0.060, oj = 54 ± 165°. However, as we show in 
Section l4~3l this is at best significant at only the 88% 
confidence level. 

We also carried out the analysis described above with 
the eccentricity fixed to zero. The resulting param- 
eters are given in Table [5] under the column heading 
b 2 ,p}, e = 0". The results are discussed further 
in Section IfHfl 

Finally, we c onducted an indep e ndent m odel of the sys- 
tem based on iKipping fc Bakosl (|2010b[ ). The primary 
differences between this model and the adopted model 
are differences in the choice of parameters to vary in the 
fit: we use T/i?* as defined in IKipping fc Bakosl (|2010bfi 
rather than CI R*i b rather than b 2 , and p 2 rather than p. 
We also allowed for a linear drift in the radial velocities 7, 
and a time shift itroj in the radial velocities due to possi- 
ble additional bodies in the system on Trojan orbits with 
HAT-P-26b. We chose to include both i tro j and 7 rather 
than fixing them to zero as the value of the Bayseia n In- 
formation Criterion (BIC; e.g. IKipping et all 120101 ) was 
lower for the best-fit model including these parameters, 
than for models where one or both of these parameters 
were fixed to zero. The resulting parameters are given in 
Table[5]under the column heading "{T/i?*, 6,p 2 }." The 
parameter values from this model are consistent with 
those from the adopted model, which gives confidence 
that our results are robust to changes in the choice of 
fitting parameters. 

4. DISCUSSION 

We have presented the discovery of HAT-P-26b, a tran- 
siting Neptune-mass planet. Below we discuss the phys- 
ical properties of this planet, and compare them to the 
properties of similar planets; we comment on the possibil- 
ity that the planet has undergone significant evaporation, 
on the significance of its orbital eccentricity, and on the 
possible presence of additional bodies in the system; and 
we discuss the prospects for detailed follow-up studies. 

4.1. Physical Properties of HAT-P-26b 

Figure [5] compares HAT-P-26b to the other known 
TEPs on a mass-radius diagram. With a density of 
0.40 ±0.10 gem" 3 , HAT-P-26b is significantly less dense 
than the four other Neptune-size planets with well mea- 
sured masses and radii (Ur anus, Neptune, GJ 436b , 
HAT-P-llb). For Kepler-4b. IKipping fc Bakosl ()2010a| ) 
find a large uncertainty on the radius which results from 
significant uncertainties on the eccentricity and the tran- 
sit impact parameter. Kepler-4b may be comparable in 
size to GJ 436b and HAT-P-llb, or it could be even less 
dense than HAT-P-26b. 

From the theoretical models of iFortnev et al.l (|2007| ). 
HAT-P-26b has a radius that is well above the maximum 
radius of 0.3 i?j for a 0.06 Mj planet lacking a hydrogen- 
helium envelope (i.e. a planet with a 100% water-ice com- 
position). The best- fit mass and radius for HAT-P-26b 
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TABLE 5 
Orbital and planetary parameters 



Parameter Adopted Value Value Value 
{(/R*,b 2 ,p} {C/R*,b 2 , P },e = {r/R„,b,p 2 } 

Light curve parameters 

P (days) 4.234516 ±0.000015 4.234515 ± 0.000015 4 -234508±°:ou0022 

T c (BJD) a 2455304.65122 ± 0.00035 2455304.65118 ± 0.00036 2455304.6512010 00049 

T14 (days) a 0.1023 ± 0.0010 0.1025 ± 0.0010 C^st^oon 

Tu = T 34 (days) a 0.0077 ± 0.0007 0.0078 ± 0.0007 0.00724± ';o 1 oo27 

a/R* 13.06 ±0.83 13 - 44 ±o.59 13. 28^°'^ 

C/R+ 21.15 ±0.16 21.14 ±0.16 21.10±o.2i 

Rp/R* 0.0737 ± 0.0012 0.0738 ± 0.0012 07341 ±oooo93 

■ 2 no92 +0 087 oiio+ 076 035- 089 

u - uyz -0.053 u - iiu -0.059 u - u,;io -0.032 

b = acosi/R^ 0.303^0^22 - 332 lo.i23 0.00±ooo 

1 (deg) 88-6lo;| 88.6 ± 0.5 89.i4to.f2 

Limb-darkening coefficients b 

at (linear term) 0.3862 0.3862 0.3862 

bi (quadratic term) 0.2576 0.2576 0.2576 

RV parameters 

K (ms" 1 ) 8.5±1.0 8.3±1.0 7.6±^;| 

fc RV c 0.099 ±0.060 0.000 ±0.000 0.09to.n 

h KV c 0.027 ±0.076 0.000 ±0.000 0.028to:oeo 

e 0.124 ±0.060 0.000 ± 0.000 °- 127 tao68 

lj (deg) 54 ± 165 ± 74 ±sg 6 

Tlms-'d" 1 ) -0.0281°;°^ 

ttroj (d) d O-Olto.23 

RV jitter (ms -1 ) 1.6 2.4 

RV fit RMS (ms^ 1 ) 2.4 3.0 

Secondary eclipse parameters 

T s (BJD) 2455307.037 ±0.162 2455306.768 ± 0.000 2455307.Ollo.29 

T M 4 0.1074 ±0.0162 0.1025 ± 0.0010 O.lOSl^on 

T s ,i2 0.0082 ± 0.0067 0.0078 ± 0.0007 - 00782 ±o.ooo83 

Planetary parameters 

M p (Mj) 0.059 ±0.007 0.057 ± 0.007 0.0522lg'^3 

R P (Rj) 0.565±°;gg 0.54S±g;g* o.553i°;^ 7 

C(M P , R p ) e 0.07 0.08 0.059 

p p (gem- 3 ) 0.40 ±0.10 0.42 ±0.08 - 378 ta084 

log 9p (cgs) 2.651^0 2.67 ±0.07 2.62lj;°;°f* 

a (AU) 0.0479 ± 0.0006 0.0478 ± 0.0006 - 04780 io.ooo6i 

T cq (K) 1001 tt b 7 981 ± 29 99ll^ 

6 f 0.012 ±0.002 0.012 ±0.002 0.0110±o 0019 

F perl (lO^rgs-icm- 2 ) s ... 2.9lt 7 : !g 2.10±g:|g 2.87± ;g| 

F ap (10 8 ergs- 1 cm- 2 ) s .... 1.81 ±0.32 2.10±q^ 1.72±qH 

(F) (10 8 erg s _1 cm~ 2 ) g 2.27tJ;° 8 2.10t°*g 2.18±g^° 

a T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. It corresponds to N± r — +31. 
BJD is calculated from UTC. T14: total transit duration, time between first to last contact; T12 — T34: ingress/egress time, 
time between first and second, or third and fourth contact. 

" Values for a quadratic law, adopted from the tabulations by Claret (2004) according to the spectroscopic (SME) param- 
eters listed in Table HI 

c Lagrangian orbital parameters derived from the global modeling, and primarily determined by the RV data. 

Time-offset in the radial velocities due to companion planets in Trojan orbits. 
° Correlation coefficient between the planetary mass M p and radius R p . 

f The Safronov number is given by = ±{V CSC /V orb ) 2 = (a/ R P )(M P /M*) ("see lHansen fc Barmanll2007fl . 
g Incoming flux per unit surface area, averaged over the orbit. 
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falls just below the 4 Gyr model with a 10 Mq rocky core 
and 8 Mq gas envelope, implying that a 4 Gyr model 
with a slightly higher core mass would provide a bet- 
ter match to the mass and radius. We note that the 
isochrone-based age (9.0^4 g Gyr) and the activity-based 

age (7. 2 Gyr) for the HAT-P-26 system are somewhat 
older than 4 Gyr, so the inferred core-mass would there- 
fore be somewhat smaller. 

We also compare HAT -P-26b to the theoretical mod- 
els of lBaraffe et al.l (|2008ft which predict more significant 
infl ation due to ir r adiati on for low mass planets than do 
the lFortnev et a l. (2007) models. In this case the radius 
of HAT-P-26b is intermediate between the Z = 0.5 and 
Z = 0.9 heavy-element enrichment models. 

4.2. Evaporation 

Observations of the transiting hot Jupiters 
HD 209458b and HD 189733b in the H I Lyman-a 
line have indicated that both plan ets are evaporating 
at a rate of up to ~ 10 10 gs _1 (e. g. iVidal-Madiar et al.1 

1\2 



10 io^g-i 



120031 : iLecavelier des Etangs et al . 2010). Prompted by 
the observations for HD 209458b, several theoretical 
studies have indicated that atmospheric evaporation is 
likely to be important for close-in planets, particularly 
those with low s urface gravities, suc h as hot Neptuncs 
(see for example iLammer et al.l 120031 and the review by 
lYelle et al.ll2008h . It has even been suggested that some 
close-in Neptune-mass and smaller planets may be the 
evaporated cores of planets which ini tially had ma sses 
comparable to Saturn or Jupiter (e.g. Baraffc 2005|). In 
the case of energy-li mited escape, the ev apor ative mass- 
loss i s given bv (see lErkaev et al.l 120071 a nd lYelle et al 



2008; see also iValencia et al.l 120101 and Uackson et al 



20101 for applications to CoRoT-7b): 



GM p K fide 



(1) 



where -Fxuv is the incident flux of extreme ultraviolet 
(XUV) stellar radiation, e is the heating efficiency and 
is estimated to b e ~ 0.4 for the case of HD 209458b 
(jYelle et al.ll2008H , and if tide is a factor that accounts for 
an enhancement of the evaporation rate in the presence 
of tides, and is given by: 



A', 



tide 



i 



3 1 

2i + W 3 



(2) 



where £ = (M v /(3M ir )) 1 / 3 a /R v is the ra t io of the Roche 
radius to the planet radius. Ribas et al. (2005) find that 
for solar type stars the XUV flux at 1 AU integrated over 
the wavelength range 1 A to 1200 A is given by: 

Fxuv.iau = 29.7r- 1 - 23 ergss- 1 cm- 2 (3) 

where r is the age in Gyr. To our knowledge, a similar 
study has not been completed for K dwarfs, however long 
term X-ray observations of the 5-6 Gyr a Cen AB system 
reveal that on average the Kl dwarf star a Cen B has 
an X-ray luminosity in the 6-60 A band that is approxi- 
mately twice that of the Sun, while the G2 dwarf a Cen A 
has a luminosity that is approximately half that of the 
Sun (Avers 20Q9|). For simplicity we therefore assume 
that the total XUV luminosity of HAT-P-26 is compa- 
rable to that of the Sun (4.64 ergs s" 1 cm" 2 at 1 AU; 



iRibas et al.l l2005h . which is likely correct to within an 
order of magnitude. Assuming e = 0.4, we estimate that 
the expected present-day mass-loss rate for HAT-P-26b 
is - 3 x lO^gs" 1 = 0.17 M e Gyr" 1 . To determine the 
total mass lost by HAT-P-26b over its lifetime, we inte- 
grate equation [T] assuming an age of 4.5 Gyr, F cx t" 1 - 23 
for r > 0.1 Gyr and F = constant for r < 0.1 Gyr, ne- 
glecting tidal evolution of the orbit, and assuming that 
the radius is constant. We find that HAT-P-26b may 
have lost a significant fraction its mass (~ 30%); the ex- 
act value depends strongly on several poorly constrained 
parameters including -Fxuv and its dependence on age 
for a Kl dwarf, e, and the age of the system. 

4.3. Eccentricity 

Using the relation given bv I Adams fc Laughlinl ( |2006D , 
the expected tidal circularization time-scale for HAT-P- 
26b is ~ 1 Gyr which is much less than the age of the 
system. This time-scale is estimated assuming a large 
tidal quality factor of Qp = 10 6 , and that there are 
no additional bodies in the system exciting the eccen- 
tricity. However, because at least two of the three hot 
Neptunes have si gnificant eccentriciti es (GJ 436b has 
e = 0.14 ± 0.01, iDemorv et aT] 120071 an d HAT-P-llb 
has e = 0.198 ± 0.046. iBakos et all 120101 the eccentric- 
ity for Kepler-4b is poorly constrained. iKipping fc Bakosl 
l2010af l . we cannot conclude that the eccentricity must be 
zero on physical grounds, and therefore do not adopt a 
zero-eccentricity model for the parameter determination. 

As discussed in Section 13.11 the eccentricity of HAT-P- 
26b is poorly constrained by the RV observations, and is 
instead constrained by requiring that the star be younger 
than the age of the universe (without the age constraint 
we get e = 0.24 ±0.12, whereas including the age con- 
straint gives e = 0.124 ± 0.060). To establish the signifi- 
cance of the eccentricity measurement, we also fit a model 
with the eccentricity fixed to zero. An F-test (e.g. lLuptonl 
19931) allows us to reject the null hypothesis of zero ec- 
centricity with only 79 % confidence. Alternatively, the 
ILucv fc Sweeney! (|1971h test for the significance of an ec- 
centricity measurement gives a false alarm probability of 
~ 12% for detecting e > 0.124 with an error of 0.060, 
or 88% confidence that the orbit is eccentric. If the ec- 
centricity is fixed to zero, the required jitter to achieve 
X 2 /-ZVdof = 1 is 2.4ms" 1 , which is closer to the typi- 
cal Keck/HIRES jitter of other chromospherically quiet 
early K dwarfs than the jitter of 1.6ms" 1 that is ob- 
tained with an eccentric orbit fit. We therefore are not 
able to claim a significant eccentricity for HAT-P-26b, 
and instead may only place a 95% confidence upper limit 
of e < 0.22. For the {Y/i?*, b,p 2 } model discussed at 
the end of Section 1331 the 95% confidence upper limit 
is e < 0.32. Further RV observations, or a photometric 
detection of the occultation of the planet by its host star, 
are needed to determine if the eccentricity is nonzero. 

4.4. Additional Bodies in the System 

The {T/R+, b,p 2 } model discussed in Section |3~3"1 with 
parameters given in Table [5] includes a linear drift in 
the radial velocities, 7, as a free parameter. We find 



7 



-0.028 



+0.014 
-0.013 



ms 



. Conducting an odds ratio 
test (see IKipping et al"1l2010f ). we conclude that the drift 
is real with 96.6% confidence, making this a 2.1er detec- 
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tion. While the detection is not significant enough for 
us to be highly confident that there is at least one addi- 
tional body in the system, this suggestive result implies 
that HAT-P-26 warrants long-term RV monitoring. We 
also searched for a linear time-shift in the RVs due to 
potential Trojans. We do not detect a significant shift, 
and may exclude |itroj| < 0.50 d with 95% confidence. 
This translates to an upper limit of 84 Mq on the mass 
of a Trojan companion, which is greater than the mass of 
the planet. With the present data we are thus not able 
to place a meaningful limit on the presence of Trojan 
companions. 

4.5. Suitability for Follow-up 

HAT-P-26 has a number of features that make it an 
attractive target for potential follow-up studies. At 
V = 11.744, it is bright enough that precision spec- 
troscopic and photometric observations are feasible with 
moderate integration times. The equatorial declination 
of S = +04°03'36.0" also means that HAT-P-26 is ac- 
cessible to both Northern and Southern ground-based 
facilities. The exceptionally low jitter will facilitate fur- 
ther RV observations, which might be used to confirm 
and refine the eccentricity determination, to measure the 
Rossiter-McLaughlin effect (R-M; discussed in more de- 
tail below) and to search for additional planets in the 
system. 

A detection of the occultation of HAT-P-26b by HAT- 
P-26 with IRAC / Spitzer would provide a strong con- 
straint on k = ecosw, while the duration of the oc- 
cultation would provide a constraint on h = esinw. 
We note that the median value of the a posteriori dis- 
tribution for the time of occultation that results from 
our global fit when the eccentricity is allowed to vary 
(Section I3.3[) is 6.5 h after the expected time of occul- 
tation assuming a circular orbit. The expected depth 
of the occultation event is a challenging 0.012% and 
0.020% at 3.6 /jm and 4.5 /im respectively. Scaling from 
TrES-4, a somewhat fainter s tar at these wavelengths, 
for which iKnutson et al.l ([2009D measured occultations at 
3.6 /im and 4.5 /im using IRAC / Spitzer with precisions of 
0.011% and 0.016% respectively, one may hope to achieve 
a ~ 1.4ct and 1.6cr detection for HAT-P-26b for one event 
at each bandpass. One would need to observe 5 and 4 
occultations respectively to achieve a 3er detection. 

Recent measurements of the R-M effect for TEPs have 
revealed a substantial population of planets on orbits 
that are significantly misaligned with th e spin axes of 
their h ost stars (e.g. iTriaud et al.l |2010|) . I Winn et"al] 
l|2010al) note that misalignment appears to be more 
prevalent for planets orbiting stars with T a g > 6250 K, 
and suggest that most close-in planets migrate by planet- 
planet or planet-star scattering mechanisms, or by the 
Kozai effect, rather than disk migration, and that tidal 
dissipation in the convective surfaces of cooler stars re- 
aligns the stellar s pin axis to the orbi tal axis of the close- 
in massive planet. [Schlaufman (2010) also finds evidence 
that planets orbiting stars with M+ > 1.2M are more 
likely to be misaligned than planets orbiting cooler stars 
using a method that is independe nt of the R-M mea- 
surements. One prediction of the iWinn etaLl (|2010aD 
hypothesis is that lower mass planets orbiting cool stars 
should show a greater degree of misalignment than higher 
mass planets due to their reduced tidal influence. The 



detect i on of misalignmen t for HAT-P-llb ()Winn et al.l 
l2010bt iHirano et al.ll2010l ) is consistent with this hypoth- 
esis. Measuring the R-M effect for HAT-P-26b would 
provide an additional test. Using equation (40) from 
IWinnl (|2010f ). the expected maximum amplitude of the 
R-M effect for HAT-P-26b is ~ 9 ms" 1 , which given the 
low jitter of HAT-P-26, should be detectable at > 10a. 

By measuring the primary transit depth as a func- 
tion of wavelength it is possible to obtain a transmis- 
sion spectrum of an exoplanet's atmosphere. Such ob- 
serva tions have been made f or a handful of planets 
(e.g. iCharbonneau et a see also th e revie w by 

ISeager fc Demintd 120101) Following |Brown| (|2001[ ). the 
expected difference in transit depth between two wave- 
lengths is given approximately by: 



where H = k^T cq / g p fi is the scale height of the at- 
mosphere, g p is the planet surface gravity, /i is the 
mean molecular weight of the atmosphere, and Njj = 
ln(eri/o- 2 ) where o\ and o-i are the opacities per gram of 
material at wavelengths in a strong atomic or molecular 
line and in the nearby continuum respectively. Assum- 
ing a pure H 2 atmosphere, \l = 3.347 x 10~ 27 kg, we find 
for HAT-P-26b H = 920km, and AS = 0.0246 N H %. 
If instead we assume that th e atmosphere has the sam e 
composition as Neptune (e.g. lde Pater fc Lissauerl feOOl). 
we have \i = 4.655 x 10~ 27 kg, H = 660 km, and AS — 
0.0177Nh%. For comparison, assuming a pure H 2 atmo- 
sphere, the planet HD 209458b has AS = 0.0198iV ff %, 
while GJ 436b has AS = 0.0107A^%, HAT-P-llb has 
AS = 0.0072A^% and Kepler-4b has A<5 = 0.0038A ff %. 
Due to its low surface gravity, HAT-P-26b easily has the 
highest expected transmission spectrum signal among 
the known transiting Neptune-mass planets. While it 
is relatively faint compared to the well stud ied planets 
HP 20 9458b and HD 189733b, we note that iSing et all 
(2010) used the Gran Telescopio Canarias (GTC) to de- 
tect a 0.058 ± 0.016% absorption level at 7582 A due to 
Potassium in the atmosphere of XO-2b, which orbits a 
V = 11.2 early K star. Scaling from this observation, 
it should be possible to detect components in the atmo- 
sphere of HAT-P-26b with N H > 3 at the - 3cr level 
using the GTC. 

4.6. Summary 

In summary, HAT-P-26b is a low-density Neptune- 
mass planet. Its low-density relative to the other known 
Neptune-mass planets means that HAT-P-26b likely has 
a more significant hydrogen-helium gas envelope than its 
counterparts. The existence of HAT-P-26b provides em- 
pirical evidence that, like hot Jupiters, hot Neptunes 
also exhibit a wide ran ge of densities. Comparing to 
the iFortnev et ail (|2007h models, we find that HAT-P- 
26b is likely composed of a gas envelope and a heavy- 
ele ment core that are ap proximately equal in mass, while 
the iBaraffe et al.l (|2008f ) models prefer a higher heavy- 
element fraction. It is also likely that irradiation-driven 
mass-loss has played a significant role in the evolution 
of HAT-P-26b-we find that the planet may have lost 
~ 30% of its present-day mass over the course of its 
history, though this conclusion depends strongly on a 
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number of very poorly constrained parameters, partic- 
ularly the XUV flux of HAT-P-26 and its evolution in 
time. We place a 95% confidence upper limit on the ec- 
centricity of e < 0.22. If further observations detect a 
nonzero eccentricity, it would mean that at least three of 
the four known Neptune-mass TEPs have nonzero eccen- 
tricities, which may imply that the tidal quality factor 
is higher than expected for these planets. Observations 
of the planetary occultation event for HAT-P-26b with 
IRAC / Spitzer would greatly constrain the eccentricity, 
however the low expected depth is likely to make this a 
challenging observation. We find suggestive evidence for 
a linear drift in the RVs which is significant at the 2. la- 
level. If confirmed, this would imply the existence of at 
least one additional body in the HAT-P-26 system. With 
an expected R-M amplitude of ~ 9ms -1 and a low stel- 
lar RV jitter, HAT-P-26b is a good target to measure the 



R-M effect and thereby test the hypothesis that low-mass 
planets are more likely to be misaligned than high-mass 
planets. The low surface gravity also makes HAT-P-26b 
a good target for transmission spectroscopy. 
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